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Abstract 

Proposed models of closed timelike curves (CTCs) have been shown to enable powerful 
information-processing protocols. We examine the simulation of models of CTCs both by other 
models of CTCs and by physical systems without access to CTCs. We prove that the recently 
proposed transition probability CTCs (T-CTCs) are physically equivalent to postselection CTCs 
(P-CTCs), in the sense that one model can simulate the other with reasonable overhead. As 
a consequence, their information-processing capabilities are equivalent. We also describe a 
method for quantum computers to simulate Deutschian CTCs (but with a reasonable overhead 
only in some cases). In cases for which the overhead is reasonable, it might be possible to 
perform the simulation in a table-top experiment. This approach has the benefit of resolving 
some ambiguities associated with the equivalent circuit model of Ralph et al. Furthermore, we 
provide an explicit form for the state of the CTC system such that it is a maximum-entropy 
state, as prescribed by Deutsch. 


1 Introduction 

Simulating one physical system using another is one of the core methods that we employ to un¬ 
derstand physics. For example, the use of classical computers to simulate physical systems is 
ubiquitous nowadays, and one of the great hopes for quantum computation is that quantum com¬ 
puters could reduce the overhead needed to simulate quantum physical systems. Thus, it is of 
interest to determine how to simulate models of closed timelike curves (CTCs) in order to improve 
our understanding of these models. 

CTCs are a tantalizing possibility not ruled out by the laws of physics, and they have been the 
subject of debate for some time now, especially since since the works in [God49[ IBonSOi IGot91] . 
Beginning with the work of Deutsch |Deu91j . the quantum information community has entered 
the discussion. The common approach has been to strip away details of the underlying spacetime 
geometry in order to focus on resolving the paradoxes associated with time travel, such as the 
well known grandfather and unproven theorem paradoxes. There are now several models under 
active consideration, including the original proposal of Deutsch |Deu91j . post-selected quantum 
teleportation (attributed to a variety of sources [BSObi ISven9[ ISvelli lLMGP'*~lla[ lLMGP'*~llb , 
see also the discussion in |Joz04j ). and most recently transition probability GTGs | All 14] . In what 
follows, we refer to these models as D-GTCs, P-CTGs, and T-GTCs, respectively. 
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Treatments of CTCs in the quantum information literature generally distinguish two subsys¬ 
tems: the chronology-violating system, which passes through the CTC and returns to its own past, 
and the chronology-respecting system, which does not itself pass through the CTC but can interact 
with the chronology-violating subsystem. The CTC itself is assumed to have finite duration: it 
comes into existence at a particular time (in the local rest frame of the experiment) and ends at 
a later time. The chronology-respecting system is prepared in some initial state before the CTC 
begins, and any measurements are deferred until after the end of the CTC. (These details are not 
essential but simplify the analysis.) 

D-CTCs are based on the following idea: for a given initial state of the chronology-respecting 
system and a given unitary evolution to take place between the chronology-respecting and -violating 
systems, nature sets the state of the chronology-violating system to be a fixed point of the evolution 
to avoid the aforementioned causality paradoxes. The resulting evolution is nonlinear: the state of 
the chronology-violating system depends on the initial state of the chronology-respecting system 
as well as the unitary interaction. 

P-CTCs resolve causality paradoxes in a different way, by modeling a CTC as a noiseless 
quantum channel into the past via quantum teleportation. Normal quantum teleportation involves 
a measurement with a random outcome BBC~*~93] . which must be communicated from the sender 


to the receiver, and hence can only arrive after it has been sent. In P-CTCs, one outcome occurs 
with certainty, and the state can be received before it is sent. Nonlinear evolution arises here as 
well because the model involves a projection and a renormalization. 

Finally, T-CTCs represent a new model of CTCs proposed recently |A1114| . with the aim of 
having fewer problematic features than D-CTCs or P-CTCs while still being somewhat plausible. 
The outcome is that they resolved a technical issue with P-CTCs, in which it is unclear what 
happens if there is no overlap between the state at the output of the unitary evolution and the 
projection. In the T-CTC model, the state of the CTC is initialized to a random pure state with 
a weight that is proportional to the probability that the evolution is self-consistent. This approach 
resolves the “null projection” issue inherent with P-CTCs. 

Several researchers have argued that D-CTCs would offer great information-processing power. 
In particular, a party equipped with D-CTCs would have computational power equivalent to the 
computational complexity class PSPACE (which contains NP) |AW09j and would be able to violate 
the uncertainty principle |BHWn9] . the Holevo bound |BHWn9] . and the cherished no-cloning 
theorem |BWW13] . In fact, recent work in YAT^IS has shown that the latter violations would 


be possible using so-called open timelike curves, in which it is not even necessary for an interaction 
to take place between the systems entering and emerging from wormholes |PRM13| . Bennett et 
al. have contested these results |BLSS09] . taking a particular adversarial model of computation 
and information processing to argue that D-CTCs would offer no added benefit for these tasks. 
This debate is ongoing [CMlOt ICMP12) . 

P-CTCs are seemingly not as powerful as D-CTCs, but still offer information-processing power 
far beyond standard quantum mechanics. In particular, Lloyd et al. have proven that a P-CTC- 
equipped party would have computational power equivalent to the computational complexity class 
PP |LMGP+llb| . a complexity class conjectured to be less powerful than PSPACE, but also con¬ 
taining NP. (The proof makes use of the results in [Aar05].) P-CTCs would also enable perfect 
distinguishability of an arbitrary set of linearly independent quantum states |BW12j , allowing vio¬ 
lation of the uncertainty principle. P-CTCs are weaker than D-CTCs, in that they cannot violate 
the no-cloning theorem |BW12] and are limited to distinguishing a linearly independent set of states 
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rather than an arbitrary set of states |BHW09] . 

An intriguing aspect of the work of Lloyd et al |LMGP^lla lLMGP~*~llb is that they con¬ 
ducted an experiment which provided a physical simulation of P-CTCs. Of course, they did not 
actually send photons back in time, but instead simulated postselected teleportation, in which they 
performed ordinary quantum teleportation, but discarded 3/4 of the experimental data so that 
the remaining data was consistent with the evolution that would occur if postselected teleporta¬ 
tion were possible. The physical simulation of generic channels from the future to the past was 
then studied in |GCH12) , which related the maximum probability of successful simulation with the 
amount of information transmitted by the channel. After these works, a different group performed 
an experiment to simulate D-CTCs |RBM'*~14) . The method to do so is admittedly simple: just 
compute a fixed point for the given initial state and unitary interaction, prepare the “CTC sys¬ 
tem” in this state, apply the unitary, and discard the “CTC system.” Of course, this approach is 
only feasible for small systems, because computing a fixed point of a given quantum circuit is a 
PSPACE-complete problem [AW09] (which is why D-CTCs have such vast computational power). 

This idea of CTC simulation points to two questions: How can we simulate CTCs using table- 
top experiments? And: Can a given model of CTCs physically simulate a different model of CTCs? 
In this paper, we address these questions in two specific ways: 


1. We prove that T-CTCs are physically equivalent to P-CTCs, in the sense that one model can 
simulate the other with arbitrarily good accuracy and reasonable overhead. This resolves an 
open question from the original paper on T-CTCs |A1114| and establishes that the information¬ 
processing capabilities of T-CTCs are exactly those of P-CTCs. As a result, the computational 
power of T-CTCs is equivalent to that of the complexity class PP. 

2. We also show how quantum computers can simulate D-CTCs if many copies of the chronology- 
respecting input state are available. Our method consists of initializing a large number of 
control qubits in a superposed state and a large number of system registers in the chronology- 
respecting input state (hence the need for many copies). We then iteratively apply a controlled 
version of the unitary interaction such that the CTC system eventually relaxes to a fixed point 
of the evolution. A final application of the unitary interaction followed by a partial trace over 
the CTC system leads to a simulation of a D-CTC. The purpose of this latter contribution is 
two-fold: 1) to offer an arguably more interesting method of simulating D-CTCs other than 
solving directly for the fixed point of the circuit, and 2) to clarify some concerns about the 
equivalent circuit model of D-CTCs [RMlOl IRD12j pointed out in |A1114| . We also provide 
an explicit form for the initial state of the CTC system such that it is a maximum entropy 
fixed-point state according to Deutsch’s maximum entropy postulate. 


The rest of the paper proceeds as follows. The next section reviews the three models of CTCs: 
P-CTCs, T-CTCs, and D-CTCs. We then present the results described above in detail. We 
conclude with a brief summary and some open questions for future work. 


2 Preliminaries 

2.1 Postselected CTCs 


We briefly review the P-CTC model [BSn5( [Sven9( ISvell( lLMGP'*~lla lLMGP'*~llb] . The chronol¬ 


ogy respecting system S is initialized to a density matrix ps, and S interacts with the chronology- 
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violating system C by a unitary Use- In the P-CTC model, the system C initially holds one half of 
a maximally entangled state, with the other half being held by a fictitious system C', which is used 
to induce a noiseless quantum channel into the past. Let denote this maximally entangled 

state: 






/c 


/C" ’ 


( 2 . 1 ) 


where d is the Schmidt rank of S interacts with C by the unitary Usc: and finally the 

state is projected onto |<k)(<h|cc'/ and renormalized. The output is thus 


1 

N 


{^Icc'Usc (ps ® |‘h)(‘h|cc") u\(j\^)cc' 


1 

N 


BsPsB^, 


where Bs = Trc {Usc}^ and one can calculate the normalization factor N as 


N = Tv{4i?5Ps} . 


( 2 . 2 ) 


(2.3) 


It is unclear what should happen in the case when the projection is null. One often-stated argument 
is that this would never happen in practice, due to slight noisy deviations from ideal evolution; 
states with null projection form a set of measure zero. Or perhaps the initial state of the universe 
is chosen such that null projections never occur. But this gap does seem to be an imperfection in 
the P-CTC model. 


2.2 Transition probability CTCs 

We now review the model of CTCs known as transition probability CTCs (T-CTCs), established 
in [All 14) . The T-CTC model supposes that the chronology-respecting (CR) system S is prepared 
in a pure state \'iIj)s and then interacts with a chronology-violating (CV) system C according to 
some unitary Use- One assumption of the model is that the CTC system C is prepared in some 
unknown state |</>)c before interacting with S. The state after the interaction is 

Usc\i’)s®\(t>)c- (2.4) 

The probability that the final state of C is consistent with its initial state is 

\\{(t>\cUsc\'^)s®\4>)c\\l, (2.5) 

where the “probability of consistency” is taken to mean the probability that a measurement after 
the interaction would find C in the same state as it started in. [A1114j then argues that the 
probability of preparation—i.e., the probability that the CTC begins in the state \4>)c —should be 
proportional to this probability of self-consistency. Thus, in this way, states that are more likely 
to be self-consistent are more likely to be prepared and vice versa. In the case that system C is 
prepared in the state |0)c) the state after the interaction is projected and renormalized to 

{4>\cUsc\'^)s®\4')c ^2 6i 

\\{(t>\cUsc\i^)s®\4>)c\\2 

Finally, since the state of the CTC is unknown, the T-CTC model supposes that it is chosen 
randomly. States could be chosen uniformly according to the Haar measure, but some of these states 
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have zero probability of being self-consistent: i.e., the quantity in (2.5) could be equal to zero for 
some \4>)ci so that such choices have zero “preparation measure.” So, by the above reasoning (that 
the probability of preparation is proportional to the probability of self-consistency), the preparation 
probability measure is biased away from the Haar measure and taken to be 


d^i{4>) = 


del) \\{(j)\cUsc\'tp)s ® \(p)c\ 


f d(f) ||(</>|c^5c|V’)s ® l<(>)c 


|2 ’ 


(2.7) 


where def) represents the Haar measure. Thus, the ensemble resulting at the output of the interaction 
is given by 

(</>|cC^5c|V’)s® l<(')c 1 ^2 8) 

(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 


and the density operator corresponding to this ensemble is 

( ) = f d iWi'il’ls ® \^){<P\c) Usc\^)c 

= {4>\cUsc ®\(l}){4’\c)ulfj\(j})c 

= ^ [(</>lct^sc|</>)c] dV') (V'Is) {(l)\cUlc\(l^)c 

with 

z = [ dej) {i}\s {4>\cUsc\(l^)c [('/>|c^5c|</>)c] IV')s- 


iPf)s = I (Trc {[/sc}\^){^I^\ sTtc {ulc} + Trc {Usc (|V')(V’|5 ® Ic) C/Jc}) > (2-13) 


|A1114| has proven that the density matrix in (2.9) is equal to 

1 

z 

where z is a normalization factor. 

In the above, if the input is the result of inputting a quantum state ps which is an improper 
mixture, then the output density matrix is given by 


(Pfh = “ {^c{Usc}psT^^c{ulc} + Tre^Usc (ps ® Ic) ^sc}) > 

where, again, z is a normalization factor. 


(2.14) 


2.3 Deustchian CTCs 

The Deutsch model of CTCs [Deu91| is rather different from the two models discussed above. The 
evolution is specified by two objects: the initial density operator ps for the chronology-respecting 
system S, and the interaction unitary Use between system S and the chronology-violating system 
C. The model assumes that the systems S and C are initially in a tensor-product state, presumably 
because they do not interact until the system C emerges from the past mouth of its wormhole. So 
the initial state after the CTC system emerges but before any interaction occurs is ps ® cre, where 
ac is a density operator for the CTC system that we will soon specify. The systems S and C 
interact according to a unitary Use-, leading to the state after the interaction 

Use {ps ® Cfc) ulc- (2-15) 
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The reduced density matrix of the chronology-respecting system after the interaction is 

Trcjf^sc (P5 ® crc) • (2-16) 

Let Nu,p denote the following quantum channel (a linear completely positive trace-preserving map): 

J^u,p{^c) = {ps ® ^c) • (2-17) 

That this map is a quantum channel follows from the fact that it is a concatenation of three quantum 
channels: the first appends system S in state ps to system C; the next is a unitary interaction Use'-, 
and the last discards system S. To determine the state of the CTC system, Deutsch has postulated 
that nature initializes it to be a fixed point of Mu,p |Deu91j : that is, a state ac such that 

J^u,p{o'c) = (2-18) 

This criterion ensures that grandfather paradoxes are ruled out. It also has the rather strange 
implication that, in a chronology-violating region, nature knows the state of a quantum system in 
a different location and also is “aware of” the unitary interaction that is about to take place, and 
then initializes the state of the CTC system based on this information. (Such effects, while strange, 
are perhaps not unexpected when information can travel backwards in time.) 

Deutsch showed that every quantum channel has at least one fixed point in the form of a density 
operator |Deu91j . In fact, the following quantum channel projects onto the fixed points of Nu,p'- 

^u,pi^c) = ^li^A/'^p {ojc ), (2.19) 

where 

_ 1 ^ 

^ ( 2 . 20 ) 

i=0 

and Mij p is defined to be i consecutive applications of the channel Mu,p- From the theory of 
operator algebras, the fixed-point set of the CPTP map Mu,p forms a *-subalgebra |HJPW04] (see 
also [MosOht rWoll2j L This implies that there exists a decomposition of the CTC Hilbert space Tic 
into a direct sum of tensor products 

= 0 ® > ( 2 - 21 ) 

3 

such that any fixed-point density operator ac of Mu,p can be written with respect to this decom¬ 
position as 

= 0 ^3 ® Pv ( 2 - 22 ) 

3 

where is a probability distribution that depends on Mu,p and the particular fixed point ac, {erj} 
is a set of states that depends on Afu,p and the particular fixed point ac, and {pj} is a set of states 
that depends on the map Mu,p but is independent of any particular fixed point ac- Finally, the 
action of the fixed-point projection map on an arbitrary state ujc is given in terms of this 
decomposition by 

= 0 Trc«^. {lijuclij} ® pj, (2.23) 

3 
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where Ilj is the projection onto the subspace T-Lcl . ®T~Lcb..- The above facts will allow us to address 
several important points regarding D-CTCs. 

The equivalent circuit model from [RMIO^ |R,D12) claimed that the map limj_^oo A/2rp(wc) con¬ 
verges to a fixed point of Nu,p and that one reaches a maximum entropy fixed point by considering 
limj^oo (^c)) where vrc is the maximally mixed state. However, the existence of a fixed point 
does not in itself imply that repeated iterations of the map converge to that fixed point. It was 
pointed out in [A1114| that this need not be the case, and an explicit counterexample was given. 
We will discuss this question in Section By invoking the above facts about fixed points, and 
providing a suitable definition of convergence, we can clarify the situation and show that indeed 
there is an equivalent circuit that converges to the D-CTC evolution. 


3 P-CTCs and T-CTCs are physically equivalent 


In this section, we prove that P-CTCs and T-CTCs are physically equivalent to each other. That 
is, each model can simulate the other up to an arbitrary accuracy, and with reasonable overhead. 

3.1 P-CTCs can simulate T-CTCs 

First we show that P-CTCs can simulate T-CTCs, using the evolution equation for T-CTCs from 


(2.14). Our goal is to exploit a P-CTC such that the final state is equal to the state in (2.14). 


To this end, we first prepare three chronology-respecting quantum registers 5, C, and C in the 
following state: 

/95 <8) [P^CC'+ (1 - P) TTC ® ttc'] , (3.1) 

where ^cc is the maximally entangled state, vr is the maximally mixed state, and p G (0,1). We 
then apply the unitary Use to registers S and C, leading to the state 


pUsc {ps ® ^cc) Ul(j -h (1 - p) Use {ps (8) TTc ® TTc) Ulfj. 


t 


(3.2) 


Now suppose that we have a chronology-violating system consisting of a single-qubit register D 
that passes through a P-CTC. We then apply the following controlled unitary: 


^ee' ® Id + {lee' — ^ee') ® ^d, 


(3.3) 


where Pauli bit-flip operator. The effect of this unitary together with the P-CTC is to retain 

the projection onto ^ee' and cancel the projection onto its orthogonal complement lee' — ^ee'- 
So after the P-CTC the final state is 


1 

zo 


p{^\ee'Use {ps <8 ^ee') \^)ce' + (1 “P) {^\ee' ^se {ps ‘Snre ^ vrc") Ul^\^)ce' 


t 


(8)4>cc' 


1 

zo 


^Tve {Use} PsTvc {ulc} + {Use {ps ® le) 


^cc', (3.4) 


_\Cf |C|^ 

where zq is a normalization factor and \C\ denotes the dimension of system C. By choosing 


p = 1/ (ICI -|- 1), we can ensure that the evolution is equivalent to that given in (2.14), which 
proves that P-CTCs can simulate T-CTCs. 
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3.2 T-CTCs can simulate P-CTCs 


We now show how to simulate a P-CTC with arbitrarily good accuracy using a T-CTC. To improve 
the accuracy requires a larger CTC system, but the scaling is reasonable. (We will show that an 
exponential improvement in accuracy requires only a linearly increasing number of qubits.) 

For the simulation, we require n single-qubit systems Ci - ■ -Cn for the CTC systems. To simplify 
the description, let |<?i>^)_R 5 be a purihcation of ps, so that TTji{\(pP){(j)f^\jis} = ps- We begin by 
preparing the state 

UscW)rs\^)cc']- (3.5) 

Next, we perform the isometry 

|0)a 0 |<l>)(<h|c'C' + |1 )a ® [Ice — ) (3.6) 

so that the resulting state is 

|0)a 0 + |1)a ® [Icc - l‘^)(‘^lcc'] Use [\4 ‘^)rS • (^•’^) 


Setting 

Po = \\\^){<^\ceUsc [\4>^)rs I^)cc']|l 2 > (3-8) 

Pi = \\[lcc' - |^)(^|cc'] Use [|</>0r 5 I^)cc']ll 2 > (3-9) 

we can rewrite the state in (3.7) as 

\p)arscc' = \/^|0)a ® \'>Po)rscc' + \/^|1)a ® \iPi)rsCC', (3.10) 

where \'ipo)RSCC' and \iPi)rscc' normalized states defined by 

\'^o)rscc' = Po ^^‘^\^){^CC'UsC [|</'^)rs|^)cc'] , (3-11) 

Hi)rscc' = Pi [Icc - l^)(^lcc'] Use • (3-12) 

Observe that the normalization factor N from (|2.3[) is equal to po- 


We finally perform the following unitary interaction between the ancilla system A and the n 
CTC systems Ci • • • 


yACi-Cn = |0)(0U ® Ici® ■ ■ ■ ® Ic„ + |1)(1U ® Xc^ ® ■■■ ® Xc„. (3.13) 


We then use (2.14) to calculate the final state: 


OC TlCi-C„{yACi-C„}\p){p\ARSCc"I^C-,-cAyACi...cJ 

+ TrCi-C„{CACi-C„(|7^)(<7’URS'CC' C) -fCi-Cn)l4ci-C„}- 
Since Trci...Cn{VACi-C„} = 2"-|0)(0|a, it follows thalQ 

TrCi-C„{VACi-C„}|¥^)(75|AR5CC'Tl'Ci...C„{PACi-C„} = 2 ^"po|0)(0|a ® |V’o)(V’o|R5Ccq (3.15) 

^We note that ( |3.15[ ) demonstrates that the vector y^|0)^ ® \'>Po)rscc' and the operator Trci•■. c„{VaCi--C n} 
saturate the bound given in IA1114I Eq. (10)]. That is, to saturate (MU Eq. (10)], we can therein set 
P = Trci---c„{VACi---c„}, IV’) = -Cw |0)/i ® IV’o)hscc', the chronology-respecting systems to be ARSCC', and 
the chronology-violating systems to be Ci • • • C„. We thank John-Mark Allen for pointing this out to us. 











TrCi-C„{VACi-C„(|</5)(‘/^URS'CC' ® ICx-Cr^yAC^...cJ = 

2"'Po|0)(OU ® |V'o)('*/’o|_R5CC" + 2"'pi|l)(lU ® I^i)('*/’i|r5CC'- (3.16) 

Thus, after normalization, ( 3.14[ ) becomes 

(2’^ + l)po 


(2" + l)po +pi 


A ® \i’o){i^o\RSCC' + 


Pi 


(2"- + l)po +Pi 


|1)(1U ® |V’i)('0i|rs'CC'- (3.17) 


With only a linear increase in the number of CTC qubits re, we get an exponential suppression of 
the undesired second term, and achieve an arbitrarily good approximation to the desired P-CTC 


evolution given in (2.2). 


4 Simulating D-CTCs 

We now discuss a few methods for quantum computers to simulate D-CTCs. The approach outlined 
here has a guaranteed fixed-point convergence and thus can be viewed as an alternative to the 
equivalent circuit model |RM10l IRD12| . We should clarify at the outset that the simulation is such 
that any third party would not be able to distinguish the output of a D-CTC-assisted circuit from 
our simulation. However, the proposed simulation need not be efficient: it could take a long time due 
to the need for convergence to a fixed point of a quantum circuit, the computation of which is known 
to be a PSPACE-complete problem [AWOQj . Note that this simulation method does not involve 
classically solving the Deutsch consistency condition as in |RBM^ . Thus, the simulation we 
propose here would be more interesting in our opinion, though far more experimentally challenging. 


Recall from Section 2.3 that the initial state of the chronology-respecting system is ps-, and the 
unitary interaction is Use- To simulate a D-CTC, it is clear that one could simulate the fixed-point 


projection of (2.20) where Ai is a large positive integer. To this end, we prepare a control register 
Fi ■ ■ ■ FA! oi N qubits in a uniform superposition of Dicke states [Dic54| of the following form: 


1 


N 




D 


N 




where 


\Dk)^ 




(4.1) 


(4.2) 


and the second sum is over all combinations of |0)®^ ® ^denoted by ® |1)®'^“^}). 

Next, we initialize N registers to A^ copies of the state ps- This requirement that we are 

able to prepare multiple copies of ps is necessary to simulate a D-CTC with an ordinary quantum 
computer—otherwise simulation is impossible. Let Ci denote the first CTC system—initialized in 
some state wci • The simulation proceeds by applying the following controlled unitary to registers 
Fi, Si, and Ci: 

| 0 )( 0 |d ® hiCi + | 1 )( 1 |fi < 8 ) Us^Ci- ( 4 . 3 ) 

Let C 2 denote the output CTC system (in general, we let Cj+i denote the CTC output of the fth 
controlled unitary). The rest of the simulation proceeds as follows: 


1. Set i = 2. 
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- output 



Figure 1: This figure depicts a D-CTC simulation for which N = 3. 


2. Apply the following controlled unitary to systems Fi, Si, and Ci'. 

|0)(0|i7’, (g) Isid + |l)(l|Fi ® Us.Ci- 


(4.4) 


3. Set i ■.= i + 1. li i < N, go to step 2. Otherwise, output system Sn (while discarding all 
other systems). 


Figure [^depicts the simulation. N needs to be large enough to ensure convergence. In practice, 
for many choices of ps and Usc^ ^ will be prohibitively large. 

The result of the procedure is the following evolution: 



Trc |c7sc (yPs ^sc\ ’ 


(4.5) 

where 

u:g ^ (u:c) . 


(4.6) 

The idea is that one could choose N large enough such that 




Trc jt^SC (^PS \usc {ps ® (Jc) 

< e 

1 

(4.7) 

for a given accuracy e > 0, where ac = M^piujc)- 

One could employ other states of the control register, such as the following state: 


1 

y/N 

(|000 • • • 0) -P |100 • • • 0) -k |110 • • • 0) -P |111 • • • 0) -P • • • -P |111 • • • 1)), 

(4.8) 


and the outcome would be the same as the above. 

We could also employ a different “iterate” for the fixed-point approximation, which might be 
easier to implement experimentally while having a faster rate of convergence. This other approach 
(called Krasnoselskij iteration [Ber07| i amounts to initializing the control register Ti • • • T/v as 
follows: 


N 



\/7|0)f, + v"! - 'y\l)Fi , 


(4.9) 


2=1 
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where 7 G (0,1) is a constant. The rest of the simulation proceeds as in the steps given above, 
leading to the following fixed-point approximation for the state of the CTC system: 


^u,pi^c) = o • • • o {ujc), 


(4.10) 


N times 


where 


= l^c + (1 - l)J^u,p{^c)- (4.11) 

By taking the limit as N ^ 00 , we find that for all 7 G (0,1) (see, e.g., [BerOTI Theorem 3.2]) 




(4.12) 


That is, the limiting map converges to the projection onto the fixed points of Mu,p- Thus, this 
procedure works equally well in simulating D-CTCs and has the advantage that the control register 
might be easier to prepare, as it is a tensor-product state. 


4.1 Maximum entropy CTC state 


The last aspect that we address is Deutsch’s maximum entropy postulate. That is: what initial 
state of the register Ci in the above simulation will guarantee that the limiting state of the CTC 
is a maximum entropy state? An initial guess is that the initial state should be maximally mixed, 
but this turns out to be correct only in some cases. 

We now show how to choose the initial state to make the state of the CTC a maximum entropy 


state. We make use of (2.23) and the formula for quantum entropy H{p) = — Trlplogp} for a 


00 

[/,P’ 


density operator p with the logarithm taken to have base two. When we input a state ojq to A/) 
the output is 

0 Trcp {nj-wcBj} ® pj. (4.13) 


Let 


qU) 


so that we can rewrite (4.13|) as 


q{j) = Trlnj-wcBj}, 


0 qU) ® Pj- 


(4.14) 


(4.15) 


The quantum entropy of the above state is 


+ '^q{j) [Hiuj) + H{pj)] , (4.16) 

j 

where H{{q{j)}) = — qU) logQ'(j) is the Shannon entropy of the distribution q. It is clear that 
we should pick ujj = nj (the maximally mixed state) because Hiujj) < H{'Kj) = logd^^-. So our 
goal is to optimize the following objective function with respect to the distribution q{j)'. 

H{{q{j)}) + ^q{j) [logdL,. + H{pj)\ . (4.17) 

j 
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We employ the method of Lagrange multipliers. Our Lagrangian is 


H{{q{j)}) + '^q{j) [logdL^ + H{pj)] + A 
j 

where A G M is a Lagrange multiplier. The derivative with respect to q{j) is 

- log q{j) - 1 + log dL^ + H(pj) + A. 

Setting this equal to zero, and solving for q{j) gives 

q{j) = 




We then choose A to normalize the (?(j)’s such that g is a probability distribution: 


(4.18) 


(4.19) 


(4.20) 


A = - log (E 2H(pj )+log dLj-l 

Thus, we can take the initial state of the CTC system to be 

= 0g(j) Pj, 

j 

with q{j) as given above, to guarantee that the state of the CTC is a maximum-entropy state. 
Also, one can easily check for this choice that 

0 {njacn^} ® Pj = ac 

j 

Presumably, for many choices of chronology-respecting states ps and interaction unitaries Use, 
it should be computationally difficult to prepare the state of the CTC system as given above. 
Otherwise, one would be able to solve problems in PSPACE efficiently. However, one can observe 
that if the second factor in each block of the direct sum is trivial, so that pj = 1, and if the size 
of each block is the same, then the maximum-entropy state of the CTC system is the maximally 
mixed state. For these special cases, the preparation of the CTC system would of course be 
computationally efficient, but the resulting D-CTC circuit then cannot accomplish any information¬ 
processing task that is not possible with an ordinary quantum channel (since, in that case, an 
ordinary quantum channel can easily simulate this D-CTC). 


(4.23) 


(4.22) 



5 Conclusion 

This paper addresses the simulation of CTC-enabled evolutions in two different ways. First, we 
prove that the recently introduced T-CTCs are physically equivalent to P-CTCs. Next, we showed 
how quantum computers with many copies of the initial input state can simulate D-CTCs, though 
not necessarily efficiently. We also addressed some ambiguities associated with the equivalent 
circuit model |HMint IRD12] , and we determined an explicit form for the state that converges to 
the maximum entropy state of a D-CTC. 
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There is an interesting and potentially practical consequence of our work here. Is it possible 
to employ the method outlined in Section and the state discrimination circuit of [BHWDQ] to 
produce a useful method for distinguishing non-orthogonal states if one has N copies available? 
To determine if the method would be competitive with existing algorithms that use a constant¬ 
sized quantum memory [BKCZ13] . one would need to determine convergence rates of the method 
outlined in Section We leave this as an open question for the interested reader. 
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